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Q . Abstract 

Q«^ ' I propose a class of D > 2 lattice SU(N) gauge theories dual to 

0^ . certain vector models endowed with the local [U{N)]^ conjugation- 

invariance and Zn gauge symmetry. In the latter models, both the 
partitition function and Wilson loop observables depend nontrivially 
^^. only on the eigenvalues of the link- variables. Therefore, the vector- 

r~| ! model facilitates a master-field representation of the large N loop- 

averages in the corresponding induced gauge system. As for the par- 
titition function, in the limit A — > oo it is reduced to the 2Dth 
power of an effective one -matrix eigenvalue -model which makes the 
C^ '_ associated phase structure accessible. In particular a simple scaling- 

condition, that ensures the proper continuum limit of the induced 
gauge theory, is proposed. We also derive a closed expression for the 
large A average of a generic nonself-intersecting Wilson loop in the 
D = 2 theory defined on an arbitrary 2d surface. 
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1 Introduction 

The theory of confinement in the D = 4 continuum SU{3) Yang-Mills 
gauge system endures as a tantalizing challenge over the last few decades. 
One possible strategy to deal with this problem was introduced in the seminal 
paper [|| of K.Wilson. The general idea is to consider a lattice reformulation 
of the continuum YM theory. In the right universality class, the lattice 
system undergoes (for certain critical value g^ of the bare coupling constant 
g^ ) a phase transition characterized by the 'divergence' of the correlation 
length. In this case the continuum YM theory (supposed to be unique in 
D = A) can be recovered as the low-energy theory at the scales much larger 
than the lattice spacing, with g^ approaching g^ along the Wilsonian RG 
trajectory. 

The employed so far D > 2 lattice YM theories are too unwieldy to 
handle analytically. Consider for example the most popular subvariety of the 
one-plaquette SU{N) lattice actions. The latter is defined by associating to 
each site x a factor 

D{D-l)/2 

E e-^««-» n XK..(t^M(x)t/.(x + /i)t/;(x + z.))t/+(x)), (1.1) 

where the multiple Fourier expansion (in terms of SU{N) irreducible repre- 
sentations R^u) involves the holonomies wrapped around D{D — l)/2 ele- 
mentary plaquettes. Besides the space-time dependence, the models like (|1.1|) 
depend nontrivially on the 0{N'^) degrees of freedom: not only on the 0{N) 
eigenvalues e'^^'^p\ j = 1,...N, but also on the nondiagonal ('angular') 
components Vtp of the relevant link- matrices Up = Up diag[e^'^^P^ f2"J', p = 
1,...,D. 

So far our analytical knowledge about systems of this type has been very 
limited, with the proper continuum limit of the D > 2 lattice gauge systems 
have been remaining beyond grasping. Besides the YM2 on a 2d surface 
§, ^ §], the available solvable examples are mainly restricted to the situa- 
tions 0, Q where the model can be transformed into some eigenvalue-theory 
of q (hermitean or unitary) matrices. The goal of the present paper is to 
take the reverse way around and induce lattice gauge theories from vector- 
field models of the eigenvalue-type. By definition, the (effective) action of 



the latter is formulated in terms of such variables that makes manifest the 
invariance under the set of a -dependent [[/(A^)]®"^ -conjugations 

s({vi/>,v^+}) = saw.}) = y: fi{Ra}) n xnAWa) , (1.2) 

{Re.} "=l 

where \I/q, a = l,..,q, is the set of arbitrary U{N) matrices. In other 
words, the largest possible conjugation-invariance (|1.2|) leads to the crucial 
simplification: the action Seff{{Wa}) ~ 0(A^^) depends nontrivially only on 
the 0{N) associated eigenvalues e*^-'*^"^ of q (unitary or hermitean) matrices 
Wa, a = l,...,g. 

To generate a consistent mapping onto the gauge models, I propose to 
start with the lattice systems defined in terms of the SU{N) link-variables 
f/p(z) and invariant under the local [U{N)]®^ conjugation-invariance 

t/,(z) ^ /i+(z) t/,(z) /i,(z) , h^EUiN) , p=l,...,D, (1.3) 

combined with the reduced gauge symmetry with respect to the center T of 
the Lie group 

U,{z)^H+{z)U,{z)H{z + p) , H{z)eT , T = Zn , (1.4) 

and, as a result of ( p^D , with the global [Zn]^ -invariance 

[T]®^ : f/,(z) ^ tpU,{z) , tp G T = Zn . (1.5) 

This class of eigenvalue-systems, being reduced to the 'one-plaquette' subva- 
riety, can be defined as following. Generalizing eq. (^[T]), one is to associate 
to each site x the factor 

E ^'' n XR,Au,i^))xR,Au.i^ + i^))xR.Au;i^ + i^))xR,Aun^)), 

(1.6) 
(where S = S{{Rfj_,^}) ) in compliance with the pattern (|T]^). What is even 
more important, owing to ( |1.3| ) the average of a (nonself-intersecting) Wilson 
loop 

Tyc(f/)=tr(^^(x)f/,(x + /i)...f/,(x-p)) = tr( J] U,,iz^) ) ^ 7^ 

{zkec} 



can be also directly expressed, as we will see, in terms of the corresponding 
eigenvalues 

<WciU)>=N'-'^c < Yl tr(f/,,(zk))> . ns) 

{zkec} 



The mapping of the model like (|1.6| ) onto the corresponding gauge system 
is to be performed in the two steps. First, we introduce the auxiliary SU{N) 
scalar field G'(z) (assigned to the lattice sites) 

f/,(z)^G+(z)f/,(z)G'(z + p) , p=l,...,D, (1.9) 

via the 'gauge-transformation' that leaves invariant both the SU{N) mea- 



sure DUp{z) and the Wilson loop observables (|1.7|). Second, one inte- 
grates over the scalar field G{z) with the Haar measure (normalized by 
J dG{z) = 1) that results in the associated effective theory with the mani- 
festly gauge-invariant action Sejf{{Up{z)}) . Indeed, after the extension (|1.9|) 
the corresponding action S{{G~^{z)Up{z)G{z + fj,)}) is invariant under the 
gauge transformations 

U,{z)=g^{z)U,{z)g{z + ^^) , G{z) = g^{z) G{z) , (1.10) 

formulated for G{z) in the unconventional way. The crucial consequence of 
the 'extended' gauge symmetry ( |1.10D is that the induced S'e//({t/p(z)}) nec- 



essarily assumes the conventional form of the multi -plaquette lattice YM 



action composed from the generic closed Wilson loops (|1 . 7|) . In other words, 
the Haar integrations over G{z) intertwine the contractions of the Up{z)- 
tensors to recollect the latter into the gauge-invariant combinations (|1.7|) . For 
example, a single yuz/ -component of the block-product in eq. ( |1.6| ) after the 
four different G{z) -integratioons ( refered to the sites of the corresponding 
elementary /iz/ -plaquette) is mapped onto 

~Xi?,.(f/.(x)f/.(x + /i)f/;(x + z/)f/+(x)) . (1.11) 

By construction, the gauge systems induced via (|1.9|) from the eigenvalue- 
models ( |1.6|) (invariant under ( |1.3|) -( pT5D ) are endowed with the same set of 
symmetries as the conventional lattice gauge theories like the Wilson's one. 
In this perspective, the local \U{N)\®^ -symmetry ( |1.3| ) can be viewed as the 



hidden' symmetry inherent (via ( |1.9| )) in the proposed family of the induced 



gauge theories. It is due to this symmetry the number of the active degrees 
of freedom per site is thus substantialy reduced 

0{{D-l)N'^)^0{DN) (1.12) 

as the action ( |1.6| ) does not depend on the nondiagonal components fip(z) 
of [/p(z) = fip(z) diag[e'^^p\'''^] fi+(z), p = 1,...,D. Nevertheless, even after 
this considerable reduction the model ( |1.6D is still not transparent enough 
to allow for an exact solution. This forces us to search for potentially solv- 
able cousins of the SU{3) eigenvalue-systems among their limAr^oo SU{N) 
counterparts keeping certain analogue of the 't Hooft coupling g'^N to be 
~ 0(A^°) . In this limit, due to ( pTSl) the eigenvalue-model (|L^) can be 
viewed as a realization of the master-field representation for the large N 
loop-averages ( |1.7|) in the associated (via (|1.9|) ) induced gauge theory. 



Indeed, in the large N limit the r2p(z) -independence of (L6) becomes 
particularly advantageous. Let the weight e"'^'^^^'"'^') in eq. ( p.. 61 ) be a generic 
function consistent with the 0{N'^) -scaling of the free energy. Owing to the 



reduciton ( p.. 121 ) of the active degrees of freedom, in the computation of the 



large A^ partitition function and loop-averages (|1.8|) one can employ the 
good old saddle-point method. The latter is to be applied either to the 
D yi{N — \) eigenvalues diag[e'''^^''^'^'>] of t/p(z) themselves (akin to [|10[) or 
to their Fourier duals (similarly to |^) - the SU{N) representations R^i,{z) 
parametrized by the corresponding sets of A^ — 1 integers Aj(/iz/|z), j = 
1, ...,N — 1 . To generate additional simplifications, we constrain the weight 
Q-s({Bf,^}) ^Q Y)e invariant under the group-product 

S{D{D -1)/2)®Z2 ; Z2: ®{p.} i?p. ^ ®{p.}i?p.. (1.13) 

It combines the group of permutations within the {/ii/}-set of the irrep- 
indices with the simultaneous Z2 -conjugation {xr{U~^) = Xr{U)) of all 
the involved irreps. The Z2 -symmetry plays, as we will see, an important 
role for the consistency of the large A^ construction. 

Altogether, the partitition function X^d of the large A^ lattice YM 
theory induced from ( |1.6|) (in a D -volume L^ ) is reproduced 



hm Xld = hm {Xrf , (1.14) 



through the partitition function Xr of the associated D{D — l)/2-plaquette 
SU{N) model with the reduced space-time dependence 






:i.i5) 



which is invariant under [Zn]"^ symmetry 

[TfD . jj^ ^ ^^^^ ^ t, G T = Zn , (1.16) 



and formulated in terms of the eigenvalues of D different link- variables U, 



p 



instead of the DL original eigenvalues of Up{z) . It is instructive to rewrite 
first the large N partitition function Xr as the specific generating functional 



{Rp\ {Rfii,} 



of the Z) -products of the generalized Littlewood- Richardson (GLR) coeffi- 
cients of the (-D — l)th order 

^ZiL} = JdU^''''''^ XR.K) [®^^ixRAUp)] e Z>o (1.18) 

assuming nonnegative integer values. 

As we will see, the invariance of e"'^*^''^^'"'^-' under (|1.13|) is sufficient for 
the reduction of the large N representation ( |1.17| ) to the 2Dth power 



limX,= lim[ V " e-s(R\D)^2D 

of an effective one -matrix SU{N) model formulated in terms of the {A} -set 
of A^ — 1 integers Xj, j = 1, ...,N — 1, which parametrize canonically the 
SU{N) irreps R = R{{X}) . The sum J2r" , being the large A^ 'image' of the 
D summations over {Rp} in eq. ( |1.17| ), is dynamically constrained by the 
following condition. Let n{R) stand for the number of boxes in the SU{N) 
Young tableau Y^^J-. associated to R. Then, in ( |1.19|) both n{R) and n{R) 
must be nonnegative multiples of {D — 1) . Next, taking advantage of the 
freedom to choose S{{Rf^u}) in eq. (|1.17| ), the Z2 -invariant effective action 



S{R\D) = S{R\D) in ( |1.19| ) can be judiciously selected in any required 
form consistent with —ln[Xr] ~ 0{N'^) . In what follows, our attention will 



be restricted to the simplest solvable class of the SU{N) or U{N) models 
with S{R\D) being defined as 



271 N 



-E E 9.(M)nE(A.-Vm (1.20) 

where g>0. In the SU{N) case the set {5'f({p})} is supposed to be invariant 
under the translations A, — >■ Aj + ?7i to match with the standard SU{N) 
condition Xjy = . As for the sum J2f i^ the exponent, it runs over the 
irreps f = r{{p}) G Y2n of the even symmetric group S{2n) (labelled by the 
partititions {p} = [lP^2P\..2nP^"] of 2n : Ellikpk = 2n) with n < Mq e 
Z>i. 

Finally, the organization of the paper is as following. In Section 2 we 
present the details of the large A^ reduction (|1.14|) representing Xr in the 



form of the GLR generating functional ( |1.17| ). In Section 3, the transforma- 
tion of the large A^ GLR functional (|1.17| ) into the 1-matrix representation 



(|1.19|) is explicitly performed, and the large A^ scaling properties of the se- 



lected 1-matrix family ( |1.19| ) are formulated. The continuum limit (CL) in 
the gauge theories induced from the models like ( |l.ti| ) is analysed in Section 
4. We propose a simple criterion for the reduced 1-matrix system ( |1.2CI| ) 
to ensure the localization {Up -^ 1} (modulo ( |1.16|) ) in the infinitesimal 



vicinity, scaling as 0{N^) , of the group- unity 1 . In turn, it predetermines 
that in the associated induced gauge theory the link-variables are as well 
localized {f/p(z) -^ 1} (modulo (|1.5|) and the gauge symmetry) which is 
tantamount to the proper CL. The phenomenon of the auxiliary 'continuum' 
limits (which accompany the conventional CL) in the effective system ( |1.19| ) 
is also discussed in connection with the large A^ phase transitions (PT). 

In Section 5 we derive a closed expressions for the large A^ average of a 
generic nonself-intersecting Wilson loop in the two-dimensional eigenvalue- 
system ( |1.6| ) on an arbitrary 2d discretized closed surface. Building on this 
expression, we associate to a given eigenvalue-system (|1.6| ) the correspond- 
ing gauge theory ( [1.1|) in such a way that both of them are supposed to 
have (in D = 2) the same CL. The general pattern of averages (including 
self-intersecting loops and irreducible multi-loop correlators) is discussed in 
Section 6. Some peculiarities of the large A^ pattern of the (multi)loop av- 
erages ( |1.7[ ) are revealed and interpreted. In the last section we make our 



conclusions and discuss possible directions of further research. The Appen- 
dices contain technical details employed in the main text. 



2 Large A^ reduction of the GLR Functional. 



Let us proceed with the derivation of the reduced representation (|1.14| ) of the 



large A^ partitition function X^d of the proposed eigenvalue-model ( |1.6| 



We will demonstrate also that reduced actions like ( 1.15 ), in fact, reproduce 
partitition function of a larger family of eigenvalue-systems including the 
'multi-plaquette' generalizations of ( |1.6| ). 

To begin with, one notes that in the theory ( |1.6|) , X^^d can be readily 
represented in the form which associates to each site z of the base-lattice 
the properly weighted sum of the blocks 



K=i / dU, ®^^; [ XBM^-u){U^iz)) Xi?,.(z)(f/p(z)) ] (2.1) 



composed of the characters which in ( |1.6| ) are refered to the plaquettes sharing 
the D p -links based at z. Employing the standard fusion- rules of the Lie 
group characters 

®Li XR^iV) = J: xrAV) 4'!|{i?.} > (2.2) 

fc=l 



the combination (|2.1| ) can be rewritten in terms of the GLR coefficients 



sD t{D-1) t{D-1) (c.n\ 

^p=l ^Rp(x)\{Rp,{x^v)} ^ Rp{x)\{Rp,(x)} ■ \'^-^) 



Consider the representation of Xj^d in terms of the GLR blocks ( |2.3|) and let 
the number of sites L^ « N while ~ln[XiD] ~ 0[N'^) . The remaining 
summations over Rfj_iy{z),Rp{z) , being parametrized by the 0{L^N) inte- 
gers, in the large A^ limit can be evaluated with the help of the saddle-point 
(SP) approach (in a manner akin to [Q; see also Section 4). Since the SP 
irreps are supposed to be space-time independent 

4°.Hz) = 4°j , RfKz) = Rf^, (2.4) 



the SP equations in the L^ -lattice theory ( |1.6| ) are reduced to the SP equa- 
tions derived from the functional ( |1.17|) associated to the reduced model 
( p..l5| ). In relating (|1.17|) and ( |1.15| ) we have used also the (7 = 2 case of 



the) identity 

^JiiSo/ n E e-^-«^''» = jhn^[ / n E e-^''"^""/^ ]^ , (2.5) 



a=l Ra a=l Ra 



valid provided that 7 > , the saddle-point values of both e ''" f '' and 

e~ ''^ '' '''"' are unique and positive, while —ln[Xp] ~ 0{N'^) . Indeed, in 
this case the equivalence of the saddle-point equations, associated to both 
sides of ( pl5|) , is sufficient to prove this identity. 

It is instructive to rederive the [X^]^ ^->- X^d correspondence (|1.14|) 
applying the SP method directly to the eigenvalues diag[e^'^^^^^^ of Up{z) . 
The effective action Seff{{ujj{p\z)}) ~ 0{N'^) , where by definition 

Xlo=JI[ ^^^^ expi-S^MuM^)})] , (2.6) 

is composed of the two parts 5*6// = Sq + Sj . The 'bare' part So{{uj}) is 
associated to the factor ( |1.6| ) rewritten as e"'^"^'''^" . The remaining contri- 
bution Sj = —Y,p,z^^[^{{^p{p\'^)})] is due to the Jacobian induced by the 
change of the variables 

N 



fdU= fdn f[ f "^A{{ijp})= fdnj{T)dT , (2.7) 



N 



A(H}) = 5[2-](^^,)ni25m(^^^)|2 , T = diagY% (2.8) 

fc=l i<i 

introducing, [7 = ^2 d%ag\e''''^\ ^'^ , the 'angular' (nondiagonal) fi- and 'radial' 
(diagonal) diag\e^'^\ -components. According to (|2.8|) and the Weyl character 
formula |0] 

X«({A})(f/) = rfetfc,,(e^^^--0/^e4,,(e^(^-^H), (2-9) 



the action S'e//({ct;j(p|z)}) (in addition to the symmetries (|1.3|) -(pr5|)) is in 



variant under the SiN) Weyl group of the the permutations 

a;,-(p|z)^^,(,)(p|z) , a^S{N). (2.10) 



In what follows, the latter symmetry is supposed to be 'fixed' by the Faddeev- 
Popov condition LUj{p\z) > a;]_|_\(p|z) . Having imposed this constraint, 

there still remains the whole orbit of the SP solutions Uj (p|z) generated 
apparently by the finite A^ symmetries ( |1.4D , (|1.5|) . 

The subtlety is that the leading 0{N'^) order of Seff{{ujj (p|z)}) (re- 
stricted to the to the SP values {u;'-'^-'} ) is invariant under the extended local 
[Zn]®"^ symmetry 

[Zn]®"" : f/p(z) ^ tp(z)f/,(z) , t,{z) G Zn , (2.11) 

which is broken by the subleading orders of Sejf down to (|1.4|),(pr5|). We 
refer to Appendix A for the details and now simply make use of ( |2.11 



In particular, it implies that (having fixed ( |2.10| )) the residual SU{N) SP 
orbit is supposed to be generated from a unique and space-time independent 
solution uj: (p|z) by the set of transformations ( p.llD which is larger than 
(|1.4|), (|1.5|) . Employing for ( p.6| ) the analogue of the identity ( p.5|) (readily 
modified to account for the degeneracy (|2.11|) ), one arrives at the required 
correspondence (|1.14|) between [Xr]^ and X^^d . 



As a side remark, in the computation of a large N average of a nonself- 
intersecting Wilson loop ( |L7|) the subleading orders of Seff{{ujj{p\z)}) must 
be necessarily included. It will be clear from Sections 5 and 6 (see also 
Appendix A). Consequently, as far as the limAr^oo < Wc{U) > averages 
is concerned, the extended symmetry ( |2.11| ) is not observable directly (al- 



though leading to certain peculiarities discussed in Section 6). 

Now we are ready to formulate the following correspondence that will be 
used in Section 4 for the analysis of the continuum limit. The latter limit 
takes place if in the (induced) gauge theory the link- variables are localized in 
a vicinity, scaling as 0{N^) , of the group-unity 1 : Up{z) — > 1 modulo ( |1.5| ) 
and the SU{N) gauge symmetry. Let in the associated reduced system (|1.15|) 
the coupling constants of Sr{{Up}) are adjusted to provide with similar 
localization f/p — > 1 (modulo ( |1.16| )). The previous SP analysis implies that 
this constraint on Sr{{Up}) ensures that in the corresponding eigenvalue- 
system ( |1.6| ) the link- variables Up{z) (entering the Wilson loop observables) 
are as well localized 

Up^i mod (1.16) ^ f/p(z) ^ i mod{(|r3),(|ri)}, Vp. (2.12) 

10 



In turn, the pattern of the mapping ( |1.9| ) guarantees that in the associated 
induced gauge theory the condition ( p.l2|) results in the required locahzation 



of f/p(z) 



Finally, let us note that the SP analysis (apphed to diag[e^'^^P^'^'>] ) re- 
veals the following direction to generalize the 1-plaquette family ( |1.6|) of 
the eigenvalue-models with the partitition function reproduced via the re- 
duced system ( p..l5|) . First, with the help of the Frobenius formula 



one is to expand the characters entering ( |1.6| ) in terms of the trace prod- 
ucts YYl=i\tr{U^{z,))Y''' . Consider a set of such products associated to the 
plaquette-factor based at a particular site x . The idea is to change the ar- 
guments in each tr(f/'^(x)) separately trading x for generic sites. The only 
natural restriction for such interchange is that the new set of the trace prod- 



ucts remains invariant under the [Zn] gauge symmetry ( |1.4| ). In this way, 
one constructs an eigenvalue-counterpart of the gauge theories with the action 
depending on the multi-plaquette Wilson loops. Provided the interchanges 
are performed in the same way for all lattice reference-sites z , the SP values 
of diag[e^'^^P^^'>] are supposed to be translationally invariant (modulo ( |2.1CI| ) 
and (p.ll|) ). Consequently, by construction of the interchange the partitition 
function of the deformed eigenvalue-system is still reproduced in the large 



N limit by the same reduced model (|1.15| ). Actually, there are many other 



ways to generalize ( |1.6|) preserving the structure of the associated reduced 
system. 

3 The effective A^ ^ oo 1-matrix theory. 

Further reduction of the GLR functional ( |1.17|) to the '1-matrix' represen- 
tation ( |1.19| ) is built on the localization of the large N summations over 
{R{p)]®{R{liv)} on the solution {i?(°)(p)}(g){i?(°)(/iz/)} of the correspond- 
ing saddle-point equations. We defer the discussion of these equations untill 
the next section, and now simply assert the expected properties of the so- 
lution in the case when the constraints ( 1.1^ ) are additionally imposed. To 



be more specific, we select the option when the effective 1-matrix system 
in eq. ( |1.19| ) is reduced to the family ( |1.20| ). Altogether, the saddle-point 

i?(°)(p) = i?(°) = i?(°) , Vp ; R^'\H = i?f = i?f , V/iz/, (3.1) 

11 



is supposed to be unique , {p}®{/uz^} independent respectively, and self dual. 
As a result, the generating functional (|1.17|) is equivalent in the large N limit 
to the reduced system obtained after the identification 

Riixu) ^R,e Fif ) , Rip) ^Re Fi^) , (3.2) 

with the remaining summations over R, R2 being localized on the same 
saddle-point values (|3.1|). 

Since the reduced action ( |1.15| ) is composed of the blocks containing equal 
amounts of Up and Up SU{N) factors, Sr{{Up}) is invariant under the 
larger set of transformations (|1.16|) where Zn is extended to U{1) . As a 
result, the partitition function Xr is invariant under the substitution of the 
SU{N) link-variables by the U{N) = [SU{N) ® U{l)]/Z^ ones 

U^UiN)^^UiN) ^ vp. (3.3) 

Complementary, the involved SU{N) irreps R^ can be viewed as belonging 
to the (anti)chiral subset of the U{N) ones. Owing to (|3.3| ), the sum in 
( |1.171 ) over SU{N) irreps R is effectively constrained by the Z2 -invariant 
pair of the U{N) conditions (nontrivial in D > 2) 

^mr-'^^ ^ <R)=<R2){D-l) , n{R)=n{R2){D-l), (3.4) 

predetermined by the pattern of U{N) GLR coefficients corresponding to 



1.18| ). In eq. (|3.4|) , the integers n{R^),n{R^) E Z>o denote the number of 



boxes in the Z2 -invariant pair of the Young tableaus corresponding to 

N N 

niRUm = Y.^^ = E(^^ -N + t). (3.5) 

Recall also that the irreps of U{N) are labelled by a set of N integers Aj 
(constrained by Yl!i=i^ K = ^n niod A^ ) 

{A^(^)} = {Ai + A^ > ... > A^_i + Xn> Xn} e [Z®^/S(Ar)] (3.6) 

generated from the SU{N) -sets oi N — l nonnegative integers {A'^^*-^-*} = 
{Ai > A2 > .. > Aat-i > 0} by the extra integer number Aat > or Aat < . 
As for the Z2 -conjugation R ^^ R, it reads 

{Xi}^{-XN-^+l+P} , (3.7) 

12 



where /?^(^) = (A^ - 1) and /5^^(^) = Ai in the U{N) and SU{N) cases 
respectively. 



Next, the effective action for R, R2 , resulting from the reduction ( p.2|) , 
contains (according to ( |1.17| )) the Dth power of L ~(gl,_i . To simplify this 

expression further, one can employ the (7 = D case of the) identity ( p.5|) . 
Altogether, it defines the one -matrix representation of the large A^ family 



hm Xr = lim [ /rff/^(^) Yl e-^^^^XR{U^) [XR.{U)f'' ?"" , (3.8) 

V— »oo W— +00 J i~z ^ ' 



R,R2 



where the weight S^^\R2) = 5'({i?(/ii^)})|{i^(^j,)=ij2} is deduced from that 
of ( |1.17| ) through the 'dimensional' reduction (|3.2|). The integrated in ( p.8| ) 



expression is simply related to the associated reduced action Sr{{Up}) . Con- 
sider first the case of the separable weight S'({-R(/iz/)}) = J2{/iu} S^R^fiu)) 
which results in S^^\R2) = D{D - l)5(i?2)/2. Once the Z2 -selfduality 
( ^.1[ ) of i?2 takes place, the identity ( |2.5| ) can be applied once more (this 
time with 7 = (D - l)/4) 

E[ e-'^^^\UU) ]'^ - [ Ee-^^^^^lx«.(f^)r ]^ = e-'^, (3.9) 
R2 R2 

where e"'^'''^''^'^^-' = e~^''^'^^''^'^\ijjp=u^ stands for the original reduced action 
( |1.15| ) with coinciding arguments. As a result, eq. (|3.8|) assumes the concise 
form 



Sr{{U}) 



lim X. = [ E jdUexvi-^^^] Xr{U^) f ■ (3.10) 



Upon a reflection, the representation ( p.lO| ) remains valid (provided the 
S{D{D — l)/2) -invariance (|1.13|) of S {{R^fiu)}) ) for a generic, not nec- 



essarily separable form of the weight S{{R{j^i')}) . 

It is appropriate to remark that the consistency of a substitution like (p.9|) 
requires the positivity Xji(o){diag[e'"^ ]) = [Xj^(o){diag[e^'^ ])]+ > 0, where 

i?2 is defined by eq. ( |3.1| ) while ul is the saddle-point orbit on which the 
integral ( pl8|) over the [/-eigenvalues Ui is localized. Also, the solution u;| 
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is supposed to be unique modulo the symmetries leaving L Jo~i invariant, 

n.\K2 



i.e. the S{N) Weyl group ( |2.1U| ) combined with the center-transformations 



Next, let us introduce the Fourier-dual partner e ^(^1-^) of e ''^d ^ 

e-«<«l-)=i:e.p|-^^lii^-L.. (3.11) 

SO that e"'^'-^^^^)/^^ = Eiji e^-^^-^^l^) XrAU) . As a result, employing the 
standard orthogonality of the characters, eq. ( p.lO| ) is readily transformed 
into the required form (|1.19|) . Note also that we have introduced an implicit 



dependence of H{R\D) on D. 

Let Z2 -invariant e~^ {R2)/2D .fg^Q^Qj- }^q g^^ arbitrary function of R2 con- 
sistent with —ln[Xr] ~ 0{N^) . The central question at this step is what are 
the general constraints defining the whole e"^^^^^^ -family of weights induced 
from the e~^ (^2)/2£>_yg^j,jg^y yjg^ ( p.ll| ). Upon a reflection, the pattern of 



the U{N) GLR coefficients ensures that the only defining property of the 
^-H(R\D) _fa^]2iily is the Z2 -invariant pair of the D > 2 conditions (|3.4|) on the 



involved SU{N) irreps. In other words, given any irrep R G ^(j:)_i),„ we 

(D-l) 



expect that there exists at least one irrep R2 G YJ^'^ so that L\^^^Jd-i 7^ . 
It suggests to factorize the constraints (|3.4]) out 



-H(R\D) _ ^^S{R\D) V^ r r _ 

e ' ' ' - e 2^ dn(R),[D-l]k dn{R),[D~l]k ■ (3.12) 

fc.fcez 

Summarizing, for a fixed D a generic residual i? -valued function e"'^*^^'^-' 
(consistent with Z2 -invariance S{R\D) = S{R\D) and with the scaling 
—ln[Xr] ~ 0{N^) ) can be induced through ( p.ll|) provided the judicious 
adjustment of e~'^*°'('^2)/2D _ 

To make the representation ( |3.12| ) practical, the two periodic Kronecker 



delta-functions are introduced via ' e -regularization' of the Poisson resum- 
mation formula 

f{n,P,e) 
'ixpy— 

kez " "pez 

f{n,p,e) = 2T[[ep'^ — i up] . (3-14) 



Y.Jn,{D--i]k = lim^exp[- p'_\J V^/[^ - 1] , (3.13) 
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Altogether, it results in the 1-matrix representation of the SU{N) GLR 
functional ( p..l7| ) 



aSToo^^ = i^J^o [[D^ T. T. T. exp{-A) r (3.15) 

L J p,pGZ nGZ>o Ji^Y^ ' 

A = S{R\D) + fHR)^P^^)+Jin{Ii),P,e)_ ^3^^^^ 



Remark that if the additional D > 2 constrains (|3.4|) in ( |1.19|) were omit 



ted, the D > 2 pattern ( |3.15| ) of Xr would be essentially two-dimensional. 



Indeed, eq. (|3.15| ) would reduce to the [2D]t/i -power of the partitition func- 



tion associated (when in eq. ( |1.20| ) g = 2 ) to the continuum generalized 
2d YM on a sphere [|, |. 

In conclusion, we note that similar analysis of the U{N) reduced model 
( |1.15| ) (with the sum running over the U{N) irreps) results in the U{N) 
counterparts of eqs. ( |3.1CI| ) and ( |3.15|) . 



3.1 The large N scaling in the selected family of the 
l-matrix models. 

Let us now discuss what conditions make the large A^ limit of the se- 
lected 1-matrix variety (|1.2CI| ) well defined and consistent with the scaling 



-ln[Xr] ~ 0{N'^), \Xj\ ~ 0{N). As we will see in Section 4, the latter 
scaling is necessary to ensure the applicability of the saddle-point method 
to the summations over the integer-valued {A} -fields parametrizing relevant 
irreps. 

To begin with, for each r{{p}) G Y2n one is to perform the proper change 
of the variables 

\- = ^j/N > 9r{{p}) = br({p})N^^ , 7^({p}) = 2~2n-J2pk, (3.17) 

fc=i 

postulating that Kdp}) ~ 0(A^°) . Next, the Weyl character formula yields 
for the irrep dimension 

dimR{{\}) = l[ {\i- \)/U - i), (3.18) 

l<i<j<N 
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which foreshadows the constraint Aj > Aj+i inherent in the pattern ( p.6| ) of 
the irrep-parametrization. Combining the latter constraint with the adjust- 
ment ( 3.17| ), one concludes that the 0(A^^) -order of —ln{Xr) is necessarily 



accumulated by those of \Xj\ which are ~ 0{N) : 

{|A,| ~0(iV)} ^^ S{{\}) ^ 0{N') , (3.19) 

which altogether justifies that the characteristic values of Xj = Xj/N are 
~ 0{N^) . For definiteness, in the derivation of eq. ( p.l9|) we have supposed 
that 5'({A})|{A =(Ar_i)/2} ~ 0(A^°) when A^ assume (unadmissible) coincid- 
ing values Xj = {N — l)/2, Vj . This choice of the additive constant in the 
definition of ^dA}) matches with the introduced pattern ( |1.20[ ). 



We will need also the alternative representation of the 'measure' in ( |1.2CI| ) 
when the constrained sum {A} G [Z®^ /S{N)] over the strictly decreasing 
integers is identically (for g > ) transformed into the unconstrained sum 
over the independent integers {A} G Z®-^ 

E <=-^«^» = 4 T. ^-^"^''. (3.20) 

where S'({Ai}) = S{{Xa(i)}), cr G S{N) . It renders manifest the 'built in' 
Weyl group of the S{N) -permutations, so that the constraint Aj > Aj+i 
can be reinterpreted as the 'fixing' of the Weyl symmetry. As a result, 
the specific scaling (|3.19|) can be viewed in fact as a consequence of the 
5'(A^)-invaraince of S'({A}) augmented by the property that e"'^*-''-'^-''-' van- 
ishes (owing to ( |3.18| )) on any boundary of the Weyl chamber, i.e. for 
V{Ai = Xj, iy^j}. 



4 The Continuum Limit and Large N PTs. 

The general idea of the continuum limit in lattice theories is captured by 
the well-known intuitively transparent condition. Namely, the properly in- 
troduced correlation length must tend to infinity (in the units of the lattice 
spacing) so that the discrete space-time is 'smoothed-out', for the low-energy 
theory, into the continuum manifold. In the lattice gauge systems, the rele- 
vant coupling constant (s) should be adjusted in such a way that the Wilson 
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loop averages < Wc{U) > undergo infinitely small changes when the con- 
tour C is deformed microscopically (i.e. at the scales of the lattice cut 
off). 

Although in a general case this criterion is rather difficult to implement 
analytically, the specific structure of the lattice gauge theories suggests a 
simpler alternative. It is intuitively clear that the infinite correlation length 



is supposed to entail the localization (modulo (|1.5| ) and the gauge trans- 
formations) of the link- variables Up{z) in the infinitesimal vicinity of the 
group-unity element 1 . To be more specific, let us fix first the 'maximal tree' 
gauge putting f/p(z) = 1 on a largest possible tree (made of the links) 
which by definition does not contain nontrivial 1-cycles. Then, introducing 
the quantum fluctuations Ap^(z) = —iln[Up'' (z)] , the required localization 
can be formulated in the large N limit in the form 

_Jim^ < [Af iz)f >~ Oif) mod (1.5) , Va, 6 = 1, ..., N, (4.1) 

where g'^N = g'^{{gk})N — i> is some 0(A^°) functional (see below) of the 
relevant coupling constants {gk} that is supposed to approach zero. The 
gauge-invariant representation of ([4.1|) evidently reads as 



hm hm jl < tr[U{pl)] > -1| ~ 0{~g'N) , (4.2) 

where U{pl) = [/^(x)[/i,(x -|- //)[/+ (x -|- z/)t/+(x) stands for the holonomy 
around an elementary plaquette in an arbitrary /iz/ -plane. 

Next, the scaling ( [4.1|) can be translated into the following (large N) 
constraint on the effective action S{{Ap{z)}) for the 'gauge fields' Ap*(z) . 
Recall that the latter can be composed of not only all possible gauge invari- 
ant oprators but also necessarily contains gauge noninvariant counter-terms 
(responsible for the restoration of the Ward identities following from the man- 
ifest gauge symmetry of the lattice gauge theory). Employing the symbolic 
form (with r{{p}) G Yn standing for the S{n) irrep similarly to ( p..20| )) 

+00 n 

siA) = E E E /?(KM), w) n M{dr^A'^)r, (4.3) 

{m} n=2 rGYn k=l 

eq. ( |4.1| ) translates into the proper scaling of the defining coefficients 

/5(r, {m}) ~ 0([A^2^]^"+^("'^™>)A^2-ELiP'=) , ^{r, {m}) > 0, (4.4) 
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where Ap = 0, Vp, is supposed to be the global minimum (modulo ( |1.5| ) and 
gauge transformations) of the action ( [4.3|) . In the case when the quadratic 
gauge invariant combination is not unnaturally suppressed (i.e. 7([2^],2) = 
but /5([2^],2) 7^ 0), the formal (large A^) 'continuum' limit A ^ of the 
lagrangian L[A) yields 

J^^o.iro^^^)^^'"^^'^'^) ' ">°' (4.5) 

where c ~ 0([^^]°) is some constant. Remark that the localization ( [4.1| ) 
does not necessarily implies that in the large N WC limit g'^N —>■ the 
effective action ( [4.3| ) becomes quadratic in F^^ . The A ^ limit in eq. 
( [4.5|) may be omitted only when 7(r, {m}) > for \/r E Yn, n > 3 . 

Finally, in D = 4 the Renorm- Group analysis ensures the existence and 
uniqueness of the low-energy YM theory renormalizable in the Dyson's 
sense (which is reassured by asymptotic freedom). That is why one may ex- 
pect that D = 4 lattice gauge theories ( [4.3| ) belong to the same universality 
class with the standard low-energy action J d'^x tr[F^^{x.)]/4:g'^j-j- . 

4.1 The {Up -^ 1} localization. 

In the particular case of gauge theories induced via ( |1.9|) from the eigenvalue- 
systems ( |1.6| ), the required localization of f/p(z) is predetermined by the 
localization ( p.l2| ) of the link- variables Up in the reduced model ( 1.15 ). For 



definiteness, the effective 1-matrix system in the representation ( |1.19D of 



( |1.15| ) is supposed to be restricted to ( p..20| ). Let us demonstrate that the 



{Up — > 1} localization ( p.l2| ) holds true provided the saddle-point values 
A- = A^A^ in ( p..20| ) approach 'infinity' according to the scaling -condition 

hm hm \Xf\^0{N/[gN-^]) ^^ fxf'^l ^ 0{[gN^-'), (4.6) 

with a functional giXg^S) approaching zero. To match with the localization 
([4. ID, ^({sffc}) in the above equation should be the same as in ([4.41). 

To begin with, we sketch under what conditions on the (rescaled accord- 
ing to ( p. 17] ) coupling constants ^t{{p]) of the) effective model (|1.20| ) the 



constraint ([4. 6]) is dynamically fulfilled. As for the overall O(A^) -scaling of 



{A'^"-'} , it has been demonstrated in Section 3. This scahng is further aug- 
mented in the WC domain when all the parameters tend to zero: {br -^ 0} . 
In this case the large N hmit of the discrete U{N) model ( |3.2CI| ) coincides 
with the corresponding hermitean model in the ezg'enwa/ue-representation 
associated via substitution ( p.l7| ) 



(4.7) 
where S{{\}\{br}) = S{{N\}\{brN^^}) + 2Ei<fc^^|Ai - Afc| . One observes 
that the sums {A} G [Zn]®^ over the independent integer-valued {A}- 
fields are transformed into the corresponding integrals over the real-valued 
variables Xj ~ 0{N^) . Complementary, in the WC domain the constraints 
( p.4| ) in eq. ( |3.15| ) are irrelevant and can be omitted. 

Defering the demonstration of these facts till the next subsection, we 
employ the WC representation ( [4.7|) to clarify the mechanism of the residual 
[g'^N]~^^'^ scaling of \\\ \ and relate the latter to the localization (p.l2|) . In 
the simplest case of the U{N) action ( [L.2(]| ) with Mq = 1 , one simply sets 



^ A^- 1 

Y^ \dimR{{X})\'^exp[-gf,J2iK-^—f] , f/2 = gf,, (4.8) 

R{{X}) i=l 

where gf^ ~ 0{1/N) and fo = [2^] . Rescaling the eigenvalues Xj = 

[gN2'\~^hj and separating in ( [4.8|) the overall factor [^^A^] * , we are 
left with the g'^ N - independent theory of eigenvalues hj that altogether 
ensures ( |4.6| ). In a general situation ( |3.17| ), similar arguments hold true pro- 



vided that all br scale to zero as in (4/?). Choosing 62/= = liminf [6^] as 
the smallest &r({p}) in the subset of r G Y2k for each particular k , one is to 
equate gN2 with the liminf [(&2A:)^] selected among all k < Mq . 

Now we are ready to make contact between the scaling ([4.6|) and the 
localization {Up -^ 1} of eq. ( p.l2| ). Recall that, according to eq. (p.2|), the 



irreps R in eq. ( |3.15| ) represent the irreps {R{p)} entering the GLR fusion- 
rules (|1.17|) . The latter ensure that the scaling- condition ( [4.6|) is valid also 
for {A*^°''(yuz/)} parametrizing the SP irreps {R^'^\nu)} on which the large 
A^ sum (|1.15|) , defining the reduced model, is localized (after integration over 



{Up} ). Combining it with ( [4.7|) , one concludes that the eigenvalues ujj{p) of 
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Up = rip diag[e^^'^''^] Vf^ are localized (modulo ( |1.16| )) in the domain where 
[Sr{{Up}) - Sr{{l})] ~ 0{N'^[g'^Nf) when ~g'^N ^ Q. According to the 
Weyl character formula (p.9|), for any particular j it results in the large N 
scaling complementary to 



lim lim < uj]{p) >~ 0{fN) mod [Z^]^ , (4.9) 



Ar^oog2^^0 



where [Zn]"^ stands for ( |1.16| ). Alternatively, in terms of the quantum fluc- 



tuations A'^p = —iln[Up''] the condition ([4 .61) can be rewritten as 

lim hm < [Af]' >~ 0{f) mod [Zn]^ (4.IO) 

for any given a,b = 1, ...,N . Altogether, it justifies the required dynamical 
localization ( p.l2| ) in the reduced model ( |1.19| )/ (|1.2CI|) under the condition 



g^N -^ . 

Let us now return to the localization ( [4 .21 ) in the gauge theories induced 
from the eigenvalue-systems (|T]^) associated to the reduced models satisfying 
( 4.101 ). The pattern of the mapping (|1.9|) (together with the saddle-point 



analysis of Section 2) suggests that the functional g{{gk}) in eqs. ( [4.1| ) and 
( [4.10| ) are to be identified. To substantiate this identification, we consider the 
large A^ WC asymptotics g'^N — » of the properly normalized partitition 
function (PF) 

X^2> = / n dUp{z)exp[-S{{Up{z)}) - S{{1})] (4.11) 

for a (induced) gauge theory on a cubic lattice with L^ sites. As the con- 
struction (|4.11| ) excludes the factor e"'^*^''^^^^ irrelevant for the contribution of 
the WC perturbative series, in the limit g'^N -^ the PF X^^ yields the 
leading order of the large A^ (g'^N) -expansion. 

Consider a generic gauge theory ( [4.3| ) with the localization ( ^.2| ) im- 
posed through the pattern (|4.4| ). Then the WC asymptotics of ( [4.11|) is 
predetermined by the ^^A^ -scaling law of the measure i.e. active Haar link- 
integrations remaining after a gauge fixing. This is simply because the exclu- 
sion of the e"'^^'^^^^ -factor guarantees that the overall action [S{{Up{z)}) — 
^di})] is ~ 0(A^^[^^A^] ) when g'^N — > . As for the Haar measure scal- 
ing, consider first a cubic L^ lattice with the free boundary conditions for 
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f/p(z) . One observes that in the hmit L ^ oo the maximal tree gauge in 
the average excludes at each site one link-integration out of the total amount 
D . The latter is clear if one imposes ('almost' complete when L —^ oo) 
temporal axial gauge keeping Ud{z) = 1, Vz . Each remaining active Haar 
integration merges, in the vicinity of 1 , with the 'flat' Lebesgue measure 
Y[p=i d^ Ap{z) . Summarizing, for L -^ oo the localization ( |4.1| ) results in 
the following large N WC asymptotics 

lim lim X[f{g) = [C ^ivl](^-i)^'^° , C > , (4 12) 

where the / dU = 1 normalization of the Haar measure is used, and C is 
a model dependent constant. Finally, one can demonstrate that the above 
pattern remains valid for a finite L provided the choice of the periodic bound- 
ary conditions for Up{z) . In the specific case of the Wilsonian lattice action 
(where in eq. ( [4.4| ) 'y{n,m) = n — 2 so that ( [4.12| ) results from a gaussian 



integration) the above asymptotics was derived in [ p!3[] . 

Summarizing, to justify the identification of the functionals g{{gk}) in 
eqs. ( ^?T| ) and ( [4.10| ), one is to prove that eq. ( |4.15| ) holds true for the gauge 



theories induced from the reduced models constrained by ( [4.6|) . To begin 
with, according to the mapping ( p,.9D , the factor e"'^^''^^^'^ can be rewritten as 
the partitition function (PF) X^o of the auxiliary model. The latter is to 
be computed in terms of the plaquatte-factor (|1.6|) where the substitution 

f/p(z)^G+(z)G(z + p) (4.13) 

is performed. Consequently, the properly normalized PF of the induced gauge 
theory can be representaed as the ratio 

xi2^=XLo/x[i (4.14) 

where X^d and X^ji are the PFs (both normalized akin to ( [4.11| )) associated 
to the eigenvalue-system ( |1.6| ) and the auxiliary model defined through ( |4.13| ) 



respectively. 

As for the large A^ limit of X^d , the correspondence ( 1.141 ) allows to 
express it as L^th power of the PF Xr of the reduced model ( |1.15|) . In 
turn, the localization ([4.10|) results in the power -like asymptotics 

lim lim X,{{gr}) = [B -gX'^f^' , B>0, (4.15) 
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which is deduced similarly to ( [4. 12| ) (while J dU = 1 is presumed). Comple- 
mentary, one easily observes that the localization (|2.12|) in eigenvalue-system 
( |1.6| ) predetermines that in the auxiliary model ( [4.13D the SU{N) field G{z) 
is localized (modulo (|1.5| )) in the vicinity of 1 as well. Consequently, akin 



to ( [4.15[ ) one obtains 



hm hm X^'^id^J) = [B ^iV^]^"^^ , B>0. (4.16) 

Combining all the pieces together, we finally arrive at the required asymp- 
totics (|4.12|) with g{{gk}) being identical to the functional which enters the 
scahng (^.6|) in the associated reduced model. 



4.2 The auxiliary 'continuum' limits and A^ ^ oo PTs. 



The 1-matrix representation (|3.15| ),( p!.20[ ) of the large N PF ( |1 . 1 7| ) allows, 



for the first time, to address analytically the D > 2 pattern of the phase 
transitions (PTs) in the SU{N) gauge theories induced from the eigenvalue- 
systems (|1.6|). In particular, the asymptotics ( |4.12| ) of Xj^2 (and similar 



power-like asypmtotics of Xj^^ -^ i^r)^ ) is characteristic of the phase 
naturally refered to as the weak-coupling (WC) one. The latter is associated 
to a connected vicinity of the WC domain {br —>■ 0} of the parameters 
defining the effective 1-matrix system (|1.20|) . On the other hand, the strong- 



coupling (SC) phase (associated to a connected vicinity of {br — >■ oo} ) can 
be shown to correspond to a completely different pattern of X j^2 {{gk}) ■ To 
take the simplest example, consider D = 2 (where the constraints ( |3.4| ) are 
trivial) and choose the simplest 1-matrix model ([4.8|) with q = 2 . As it is 
known |jl2[ , the large N free energy of the latter model can be expanded in 
the SC series 

+00 M{n) 

- lim ln[Xi{b)] = N^Y.Y. /(^'^) b^""exp[-nb^/2] (4.17) 

iV ^oo 

n=l m=0 

which reproduces the correct answer in a connected vicinity of 6 = g^N — »• 

oo . 

The reason for the mismatch between the WC and SC large A^ patterns 
resides in the phase transition(s) which are ubiquitous p, |^ in the lattice 
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gauge systems. As we will see, in the case of the induced gauge theories, 
it implies that the large N continuum limit (CL) encoded by (^l3| ) - (|4.5|) is 
accompanied by a few auxiliary 'continuum' limits in the effective 1-matrix 
system (|3.15|) . The latter are characterized by the transformation of a few 
relevant discrete spaces into the corresponding continuum manifolds. 

As previously, we concentrate on the simplest case when the 1-matrix 
model in (|3.15| ) reduces to the selected family ( p.. 201 ) . The appropriate frame- 
work for our analysis is provided by the method based on the saddle-point 
equations (originally introduced in |Tl|) in a particular realization close to 
@, 1^. First of all, let us demonstrate that in the large N limit the SU{N) 
D > 2 constraints (|3.4|) in eq. ( p.l5 ) are irrelevant and can be omitted. The 



easiest way to see it is to notice first that for a self dual SU{N) irrep 

n{R)\^^R = NiX,-N + l)/2 (4.18) 

so that, choosing A^ G 2{D — l)Z>o , we satisfy ( p.4|) for V{A} . Presuming 
the 'smoothness' of the large A^ limit, one arrives at the required conclusion. 
Actually, employing the 'regularized' Poisson resummation formula ( p.l3|) . 



the irrelevance of ( ^.41 ) in the SU{N) models can be explicitly proven with 
the help of the formalism we now focus on to handle the effective 1-matrix 
system ([T^ ). 

As we demonstrated in Section 3.2, the adjustment ( |3.17] ) of the {gk}- 



scaling is sufficient to ensure that the finite A^ sum ( |1.20|) is accumulated 



by the configurations {Xj} possessing a 'smooth' large A^ limit both in the 
index j -space and in the 'base' {A} G Z®^ -space. Then, in the computa- 
tion of the leading 0{N'^) order of —ln{Xi) the large A^ sum ( p.. 20] ) can be 
substituted by the corresponding 'path-integral' 

Xi = j DXit) exp{-N^Seff\X{t)] ) (4.19) 

A(1-^) = ^ = A. , N'S^ff[X{t)]^S{{X]\ (4.20) 

over the ' continuum ' variables X{t) evidently obeying the inequality |^ 
following from ( p.6| ) 



A. - A,+i - ^ > 1 , t = {l- r/N) G [0, 1] . (4.21) 
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If the parameters {gk} are such that the solution \^^\t) of the saddle-point 
equations 

5S^jf[m] ^ ^ ^ A^^W(o) = o, (4.22) 

oX{t) 

fulfils dW{t)/dt > 1 for Vt G [0,1], the U{N) pattern of ( ^l9D reduces 



to that of the large N hermitean model associated to (|1.20|) via the corre- 



spondence ( |4.7| ). As we will prove below, this reduction occurs in the WC 
phase when the set {gr} = {K} belongs to a connected vicinity of {br -^ 0} . 
In the corresponding hermitean model, the 1/A^ quasiclassical expansion is 
organized in the background of the saddle-point solution(s) \^'^\t) . Since 
S'({A}) ~ 0{N'^) , the quantum fluctuations of the continuum A(t) -variables 
can be neglected in the large N limit 

^irn^Xi = exp{-N'S,ff[X^'\t)]) (4.23) 

and the leading order of the semiclassical and that of the 1/A^ expansions 
coincide . Remark that, in the case associated to the SU{N) system, we are 
to impose the additional 'boundary' condition ( |4.22| ). The latter matches 



with the translational invariance Aj — >■ A, -|- m which must be present in the 
(Z2 -invariant) SU{N) action ^({A^^(^)}) defined by eq. (CT) . 

One can show that the SU{N) solution of (f4.22|) , (|1.2CI|) does exists (pro- 



vided the subset {gr, r G Y2M0} has the signs consistent with the convergence 
of the {A} -series), being unique and Z2 -invariant. In the U{N) -case the 
above m -translations (of Aj ) generate the 'quasizero' mode (to be factorized 
out in the analysis of the 0(A^^) -order of — /n[Xi] ). In effect, it suffices to 
retain the (unique) Z2 -invariant U{N) solution which satisfies 

/i(o)(t) = -/i(o)(l - t) , h{t) = X{t) - 1/2 , (4.24) 

in terms of the shifted function h{t) . The explicit construction is easier to 
perform reformulating ( |1.20| ) in terms of the spectral density associated to 



/iW(t) 

r-l 
/O 



< piv) = I dt6{r^ - h^'\t)) = (rftM(°)(t))U(o)(,)=, < 1 , (4.25) 



constrained by j drip{ri) = 1. Remark that, in the case corresponding to 
the SU{N) model (|1 . 20| ) , one is to impose additionally that p^'-^^^'{r]) = 
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if V?7 < — 1/2 to match with the condition ( [4.22| ). Note also that condition 
( [4.21| ) results in the constraint pirf) G [0, 1] . 

To take the simplest example, consider the U{N) subvariety of (|1.20|) 
when the factor in the exponent e~^*^'^^ can be rewritten in terms of h{t) in 
the 'local' form 



iV2 



jdtV{h{t)) , n^)=E ^ ."^^ > (4-26) 

-' m=l ^"^ 



corresponding to p2k = 1, P2k~i = 1, 1 < ^ < Mq . In this case, the SP 
equation for p{ri) reads 

l''VM_pf^m^ . p(^,=p(_^), (4.27) 



q dr] J 1] — 

where P stands for the principle value of the integral over , and the ad- 
ditional reflection-invariance condition selects the solution associated to a 
selfdual SP irrep R^^^ = R^^^ . There exists a standard algorithm to solve 
( [4.27| ) for an arbitrary potential ( |4.26D . It can be found e.g. in the third 



reference of 0] and will not be reviewed here. We note only that, in the 
case of the simplest Mq = 1 U{N) model (f4.8|), the solution p{ri) of P3 
assumes the pattern of the Wigner semicircle law 




Piv) = -\R-V' , b^b,, (4.28) 



71^ 



which allows to reproduce the original descrete sum ( |1.20| ) for (26^/g) < 
(when p{rj) < 1 in the whole admissible domain of r]"^ < 2g/6^ ). 

Next, in the domain of {gr} where dX^^\t)/dt < 1 in a set of 'windows' 
t G Ufc [ia{k),ib{k)] , the above strategy has to be refined. To be more specific, 
consider again the simplest Mq = 1 U{N) model ([4.8|) (analysed at length 
in 0). When 6^ > 6^r = Q'7r^/2 , the solution X^^\t) is located partially 
outside the admissible domain ([4.21|) as it is clear from ( [4.28| ). To reproduce 
the original discrete model ([4.8|) with 6^ > b^^ , in ( [4.23| ) one is to substitute 



A'^°^(t) by the 'closest' boundary configuration C^^\t) . The latter is to be 
found presuming that 

2^ = 1 , te[t,M, (4-29) 
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where the [Z2] -sefduahty prescribes that — ti + 1/2 = t2 — 1/2 • The choice 
( [4.29| ) imphes simply that, for the indices j corresponding to the window 
[^1,^2] , the number Uj of boxes in the associated jth rows of YJ^^'' (con- 
strained by rii > nj-|_i ) assumes the 'boundary' values rij = -n-j+i = .. = 
riM-i = 0. After substitution of ( |4.29| ) into the effective action ( [4.20| ), one 



obtains a (modified compared with ( |4.22|) ) saddle-point equations for C^'^^t) 
in the remaining domain t G [0, ti] U [t2, 1] ■ 

Practical computations are easier to perform translating the above con- 
struction in terms of the spectral density ( [4.25| ) to be deformed p{ri) ^ 1 on 
the interval rj G [i]i,ri2] corresponding to ( [4.291) . The 'turning on' ( [4.29| ) of 
the constraint (f4.21|) (being a nonanalytical procedure) results in the third 
order phase transition p. Note also that in the Mq > 1 case of (|1.2CI|) 



there may be multicritical pattern with a few phase transitions associated to 
opening or closing of a new 'window' [ta(i),tb{i)] (see e.g. the third ref. in 

i). 

Next, the above picture of PTs allows to reinterprete the mismatch be- 
tween the large N WC and SC patterns of the partitition function X^d 
(eqs. ( |4.12| ) and ( [4.17| ) respectively) from the viewpoint of the auxiliary 



continuum limits. To begin with, according to the identification ( 4.21| ) the 



N ^ 00 limit of the discrete index-space {i = 1,...,A^} is the continuum 
'time '-manifold t . Outside the WC domain, the large A^ system ( |4.19| ) still 



retains the remnant of the discreteness (associated to the base Z®^ -space) 
which is encoded in the constraint ( [4.21| ). The subtlety is that the latter 



constraint is relevant despite the fact that each particular large N sum over 
a given j -species Xj G Z merges with the integral over the real- valued vari- 
able. Finally, let us show that in the WC limit {bj. -^ 0} the solution X^^'{t) 
doesn't violate condition ( [4.21| ). In other words, all the involved discrete 



spaces merge with the associated continuum manifolds which, in particular, 
implies the validity of the large A^ reduction ( [4. 7] ) of ( |1.2(j| ) in the whole 



A^ ^ 00 WC phase. Indeed, the substitution ( [4.7D can be justified by going 



over to the effective theory of the Xj = [9^^]^j fields. Having adjusted the 
proper scaling (|3.17| ), (^^ ) (i.e. ^^A^ -independence of the large A^ Xj -system 
in the WC limit), one obtains 

dX{t)/dt = [gN^-^dx{t)/dt ~ 0{l/[gm]) » 1 (4.30) 

which substantiates the WC equivalence of (|1.2CI|) with the associated her- 
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mitean model. 

In conclusion, it takes a minor modification of the above analysis to retain 
explicitly the D > 2 SU{N) constraints (|3.4| ) built into the effective 1- 
matrix model ( |1.19| ). Employing the representation ( |3.15| ) of the latter, one 



can demonstrate that the free energy N'^Feff{p,p) (evaluated prior to the 
sum over p,p E Z ) assumes the form 

jhn e-^^[^^//(^'^-)-^^//(°'°)l = lim e^^^ 5,,o E ^p,[i^-i]^ • (4.31) 

As a result, the localization {p = 0; p = mod [D — 1]Z} proves that the 
constraints (|3.4| ) can be safely omitted in (|1.19| ), (|3.15| ) when N -^ oo . 



Summarizing, the large A^ continuum limit (CL) in the proposed induced 
lattice gauge theories is accompanied by the two auxiliary 'CLs'. The latter, 
being predetermined (in the reduced system (|3.15| )) by the correspondence 
( p.l9| ) and the scaling ( ^.6| ), are associated to the N ^ oo transformation 
of the discrete j -index space and {A} G Z®^ space into the corresponding 
continuum manifolds. In turn, it foreshadows the irrelevance of the D > 2 
constraint ( p. 41) in the limit N ^ oo . Together with the large N reduction 
( p..l4| ), it results in the substantial simplification of the large N WC analysis. 
As for the physical interpretation of the discussed phase pattern, it can be 
understood relating the large A^ PTs to the issue of the validity of the 
WC and the SC expansions in the associated WC {br — > 0} and the SC 
{br -^ oo} domains respectively. We expect also that, among the phase 
transitions (passing from the WC to the SC limits), there is one associated 
to the condensation of the microscopic lattice monopoles which are always 
present in the lattice gauge theories. 



5 Remarks about the D=2 case. 

To begin with, consider the D = 2 large A^ average of some nonself- 
intersecting Wilson loop ( |1.7| ) for which a closed expression can be derived. 
The appropriate method is a minor modification of the combinatorial ap- 
proach due to Migdal (see also 0, ||) developed in the context of the 
ordinary 2d gauge theories. 
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Owing to the periodic boundary conditions the 2d L"^ -plane in question 
is made, topologically, into a discretized 2-tora. Let L — i> oo so that any 
finite-size loop C does not coincide with one of the two uncontractible cycles 
defining the tora. Consequently, cutting along the contour C , this tora is 
decomposed into the two disjoint connected componenets: a disk with area 
A and the complementary disk with one handle of the area L"^ — A . To ap- 
ply cutting-gluing technique, first one is to introduce (as in ^j) the so-called 
'disc-amplitudes' Zg{A, Lc\{U{zi^)}) . The latter is defined as the partitition 
function of the D = 2 eigenvalue-theory (|1.6|) on a disk of the area A with 
a given number of handles G . The boundary-contour C is supposed to con- 
sist of Lc pfc -links, each being endowed with a free boundary conditions 
introduced via the associated link-variables f/pj.(zk), k = 1, ..., Lc ■ 

A particular term of each defining factor (|1.6|) , being refered to a fiu- 
plaquette based at z , assigns to this plaquette one of the SU{N) irreps 
R^yiz) = -R(z) . The specifics in the computation of the 2d disc-amplitudes 
is that Zg{A, Lc\{U{zk)}) can be evaluated with the help of the orthonor- 
mality condition for the Lie group characters 

LrI\R. = ldV XrAV) xrAV^) = 6r,,r,, (5.1) 

corresponding to the simplest GLR coefficient of the second order while the 
higher order Lr ,cj^ i, p >3, are not involved. Given the amplitude associ- 
ated to a particular connected base-surface, eq. ( ^TTl) prescribes that R^ui^) 
is z - independent . As a result, thus introduced disk- amplitude is readily 
computed in the G - independent form 

^G(A|{f/(zk)}) = Ee~"^''^^ n Xij(f/p.(zk)) ^5 2) 

R {zkec} ^ ^ 

where {zk} stands for the set of Lc sites on the boundary C . As a cross- 
check, the D = 2 partitition function Xl2 can be easily rewritten in terms 

of ZG(A|{f/(zk)}) (where A^ = A, A2 = L^ - A) 

Xl^ = JI[ dU,,{z^) Zo(Ai|{f/+(zk)}) ZM2\{U{z^)}) (5.3) 

{zk} 



which, owing to ( |5l^ ) and ( p.5| ), results in 

r2 , ^-^ S(fl) 



tl^^ = ^i" S»"'""' = I ^'-^ r (5.4) 



H H 
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that matches with the D = 2 case of ( |1.17| ) according to ( p.l| ). Remark that 



the same expression comes out in the case of the L^ -surface of an arbitrary 
genus G . 

Next, similarly to the ordinary 2d lattice gauge theories, the Wilson loop 
average in the 2d system ( |1.6|) reads 

XRoiUc) 



< W^c" >= / n dU,,iz^)ZoiA,\{U^iz^)})ZM2\{Uiz^)}) ^ , 

•^ {zk} ^L^ 

(5.5) 
where Uc = Ilizk} ^ptl^k) , Ro is the representation of the holonomy. Eval- 
uating ( p. 51 ) one is to take advantage of the fact that ZG(A|{f/(zk)}) is 
independent of the 'angular' variables ^(zk) . Equivalently, it can be repre- 
sented as the invariance of ZG{A\{U{zk)}) with respect to the [U{N)]®^'^ 
conjugations 

f/p.(zk)-V;t(zk)f/p,(zk)l-,,(zk) , V,,iz^)eUiN), (5.6) 



(where k = 1, ..., Lq ) inherited from the symmetry ( |1.3| ) of the factor ( |1.6| ). 
Let us perform the substitution (|5^ ) as a change of the variables in eq. (^.5|) 
that leaves the measure invariant. Integrating over V^^(zk) with the Haar 
measure and employing the identity 

/ n dviz,)xn.{Y[ vHz,)u,Mv{z,)) = "^;:;.^^;^!^:^r^l (5.7) 

one trades the factor XKodljzk} Up^{z]^)) for the properly normahzed product 
of the characters (akin to ( [L.2|) ). At the same time, the rest of the integrated 
in ( |5.5| ) expression is left intact. Remark that in eq. ( p.7| ) the [dimRf)Y~^'^ 
normalization can be deduced making the particular choice of the arguments: 

[/,,(zk) = i, vfc. 

After {dV{z.p)} integrations, a generic loop average (|L^) in the eigen- 
value systems can be expressed in terms of the GRL coefficients ( p.2|) . For a 
nonself-intersecting contour C one obtains 

^ wi?.o ^_ J_ V p-S{Ri)A-S{R2){L^~A) [Lr,\R2,Ro]^'' , . 
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which complements the GLR pattern of the partitition function. Observe 
that the same GLR pattern (|1.8| ) of < W^ > immediately follows from (p.7|) 
in any D >2 . Finally, < Wq° > can be rewritten in a more conventional 
form of the average 

< W^c° >= ^E^"'''''^''^ Bl^{R\Ro,A) = < BLa{R\Ro,A) >R (5.9) 
performed with the weight e~^ ^^^^ , and we have introduced 

]Lc 



BL,iR\Ro,A) = ^e-[^(^^)-^(^)]^ j^ 



[dimRol^c-'^ 



(5.10) 



Remark that eqs. ( |5.9| ),( pTT0D remain valid in the case of a generic nonself- 



intersecting contour C on a discretized 2d surface of an arbitrary genus 
(provided C is not homotopic to any of defining uncontractible cycles). 

Building on the results of [Q (obtained in the context of the continuum 
YM2 on a 2d sphere), one can derive the integral representation for ( ^.9|) in 
terms of the associated spectral density ( [4.25] ). For simplicity, we consider 
the option when e"'^*^'^-* belongs to the U{N) subvariety (|4.26| ) of the family 
( |1.2(J| ). Also, we restrict our attention to the segment of the WC domain of 
{^^} = [br} where the spectral density is less than 1/2 for all admissible rj . 
Then, as it is demonstrated in Appendix B, the large N average < W^ > 
of a nonself-intersecting loop C in the fundamental representation assumes 
the form 

< wi >= iv-'c |,,„ ,(„) (^H!^y' , (5.11) 

where p{ri) < 1/2 is the Z2 -invariant solution of the SP eq. (|4.27| ). We 



postpone the interpretation of the unconventional A^^ ^^ -scaling of the av- 
erage ( [5.11| ) till the next section and now turn to the issue of the continuum 
limit (CL). 



5.1 Identification of the universality class. 



To study the CL of (|5.11|) directly, one is to send L'^,A to infinity keep- 



ing A/L'^ finite. Simultaneously, the coupling constant b^ = g^N (entering 
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LGD ) should approach the critical value b^ = (corresponding to the lo- 
calization (|2.12|) ). It is to be performed in accordance with such scaling-law 



which ensures a finite physical string tension a as 6^ — *■ . Let the total 
area of the closed 2d surface in question, being measured in the units of 
a^^ , is adjusted to be finite as well. Take the G = case, the 2d sphere, 
where the free energy of continuum YM2 is ~ 0{N'^) . Then thus defined 
6^ ^ limit of the (area-dependent part of) the G = amplitude (|5.11| ) is 
to be compared with its continuum counterpart. The latter can be directly 
computed with the help of the self -reproducing lattice theory [P, |^ defined 
via the plaquette-factor ( |1.1|) with 

e-^(^) = dimR exp[-f{{Ck{R)})] ■ (5.12) 

Here f{{Ck{R)}) is some function of the Casimir operators Ck{R), k = 
1,...,N, which parametrizes a particular generalized YM2 [§]. Note that 
the q = I case of the family ( p..2U| ) represents the polynomial subvariety of 
( ^.12| ) for the particular choice of U = 1 . 

The identification of f{{Ck{R)}) , corresponding to a given eigenvalue- 
system (|1.6|), would determine the universality class to which the associated 
induced lattice gauge theory belongs. Instead of addressing this question 
directly, we will find a lattice gauge theory (|1.1| ) which is supposed to have 
the same CL (in a finite physical 2-volume) as the model ( |1.6| ) in question. 
For this purpose, in the infinite lattice- volume limit L^, A — > cxd we identify 
the pairs of the large N gauge- and eigenvalue-systems for which the large 
A loop-averages are the same for generic values of the relevant coupling 
constant (s) {br} . 

To begin with, from the results of 0, 0, one readily deduces for the 
disc-amplitude associated to (|1.1| ) 



Zg{A\Uc) = Y.idimRy-^'' e-^(^)^ xr{Uc) , (5.13) 

where K{R) = F{R) -\-ln{diniR), . To make a crosscheck, one observes that 
for the Heat-Kernal action (^ (with f{{Ck{R)}) ~ C2{R) ) eq. (|57[^ ) re- 
produces the canonical expression derived in [|, || . Compared to the pattern 
( p. 21 ), the above formula contains the topological piece {diniRY~'^^ (which is 
A -independent) and depends on the eigenvalues of the single holonomy Uc ■ 
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As a result, the large A^ partitition function Xj^2 on a discretized surface 
of genus G assumes the form 

jm ^(rf^mi?)2-2G e-^(^)^' = [ ^(d^mi?)^ e"^ ]2^' , (5.14) 

R R 

where the weight K{R) is supposed to be consistent with the 0(A^^) -scaling 
of /n[X|2 ] • To be more specific, in what follows we concentrate on the case 
when e~^^^^ in eq. (|1.1|) is restricted to the family ( |1.20|) . For a preliminary 



determination of the pairing, one observes that the pattern (^.14|) precisely 
matches with that of eq. (^.4|) provided the identification 



F{R) + ln{dimR) = K{R) = S{R) , L^ ^ 00 , VG , (5.15) 

complemented by the infinite lattice- volume limit. To ensure that the free 
energy (in both ( |1.1| ) and in the associated via ( p.l5| ) eigenvalue-system (|1.6|)) 
is ~ 0{N'^) , one is to require that the pattern ( p.. 20 ) of e~^^^^ satisfies 



q > max[{l - (2 - 2G)/2L'^) ; 1] . Note also that, in the G = 1 case of the 
2-tora, the matching (|5.15|) accidentally holds true for any finite L^ as well. 

Next, consider an arbitrary nonself-intersecting loop C of the area A on 
a generic discretized 2d surface of genus G with the total area L^ . For 
simplicity, we restrict our attention to the case when, after cutting along the 
contour G , the disc of the area A has zero number of handles. Employing 
0, H, one derives 

< W^ >= -J^ Y.idimRf-''' e-^'^(^) B{R\Ro. A) , (5.I6) 



X 



(G) 
L2 R 



B(m<>, ^) = E §^ ,-iKm.yKi.M i^_,^^^_^ . (j^^^j 



Ri 



dimR 



To simplify the comparison with ( |5.9| ), in what follows Ro is selected to 
be the fundamental irrep / . According to the representation theory , it 
substantially reduces the pattern of the GLR fusion rules 

{Ln^^nj = 1 or 0} =^ {Lr,\rj)''^' = Lr^irj (5.18) 

which makes the relevant analysis particularly transparent. Combining (|5.18|) 
with the identification ( |5.15|) , one observes that the peculiar 'renormaliza- 
tion' of the perimeter-law by the factor {dim /)^~^^ = N^~^^ is the only 
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difference between ( p.9| ) and ( |5.16| ) in the limit L"^, A ^ oo . More precisely, 
denote the effective actions in the integral representation (akin to eq. ( ^.11|) ) 
for the loop-average 

< Wl >= N^'-^'^o) jdh e-^'/'/C'l^'i''^}) (5.19) 

by B).^j{h\A,{hr}) where k = 1,2 stand for the options ( |5.9| ) and ( |5.16| ) 
respectively, and Hi{L) = 2 — L, H2{L) = 1. According to the exphcit 
form ( [OOD of e-^(^) = e'^^^) , the large A limit lim^^oo Siff{h\A, {6,}) is 
/c -independent for generic values of {br} . To be more specific, consider the 
continuum limit 6^ — i> and take the Mq = 1 U{N) model ( |4.8| ). In the 
limit 6^ — > cx) the spectral density ( |4.28D tends to zero: p{ri) r^ 0{b) . As it 
is clear from Appendix B, this scaling of p{ri) implies that [dimRi/dimR] — >■ 
1 -|- 0(6^) (for admissible r]"^ < 2q/b'^ and the characteristic values of R). 
In other words, the [dimRi/dimR] -dependent mismatch between ( p.lO| ) and 
( ^.17| ) is expected to be irrelevant in the continuum limit. 

Upon a reflection, the above correspondence between (|5.9| ) and (|5.16|) 



suggests the following equivalence. The gauge system (on a sphere of a 
finite physical 2-volume) induced from ( p..6|) and the associated ordinary 
gauge theory (|1 . 1| ) in D = 2 have one and the same continuum limit (CL) 
provided the identification ( 5.15|) . To justify the general consistency of this 



statement, we observe first that once the localization (|2.12D takes place in 



the 2d system ( |1.6|) then similar localization (with the same scaling) is valid 
in ( |1.1|) associated via ( |5.15|) . The converse statement is also true. This 



correspondence between the localizations follows from the specific relation 
( ^.15| ) between the pair of the weights e~^^^^ and e~^^^'^ . To make it man- 
ifest, one simply needs to change the variables going over from Xj- to the 
Xj = {[gN^Xj/N)-Re\ds like in eq. (g]^. 

Next, the localization in the ordinary D = 2 lattice gauge theory indeed 
implies that the latter theory approaches the continuum. It is predetermined 
by the absence of the propagating degrees of freedom in the 2d gauge systems 
so that the renorm-group anomalous dimensions vanish. To see how it works, 
introduce the infinitesimal lattice spacing a and some finite parameter 
K ~ ^/a responsible for the physical mass-scale. Then, given the localization- 
pattern (|4^ ) of the lattice action ( ^73|) , the judiciously adjusted scaling 



Ap = Bpa , h = na ; lim b L = const , (5.20) 



+0 
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in the limit a ^ is supposed to directly provide with the pattern of the 
continuum gauge action associated to the corresponding low-energy theory 

of(0)- 

Finally, the equivalence of the CL in ( |1.6| ) and ( [1 . 1| ) is tantamount to the 
equivalence of the effective theories for the low-frequency modes of Bp . The 
perimeter-type Lc -dependent terms, among other things, are responsible for 
the details of the s/iort -distance attachment of the colour-electric flux to the 
external Wilson loop source. As such, they may be sensitive to the details 
of the discretization (of both the space-time and the action) that in our case 
is reflected by the mismatch in Hk{Lc) of eq. ( ^.ISI ). On the contrary, in 
the CL limit { WM (keeping A/L"^ finite), the bulk (fe^L^), (6M) -dependent 
terms in both ( ^^ and (|5.16| ) are supposed to be low-energy quatities insen- 
sitive to the lattice discretization. 

Formally, one could consider instead of a single loop average the so-called 
Creutz ratios (see e.g. [l^) of the products of the averages 

L+ = Lci + Lc2 = Lc^ + Lc4, (5.21) 



< Wi X Wi > 



where the pairs of the contours {Ci,C2} and {0^,04} are constrained to 
have the same overall perimeter L+ . Being originally designed to cancel 
exactly the short- distance perimeter-type contributions, the ratios are evi- 
dently insensitive to the difference in Hk{Lc) of eq. (|5.19|) . Altogether, the 
above general arguments are consistent with the presumable equivalence of 
the universality class to which the latter pairs of the ordinary and the induced 
2d gauge theories belong. Note also that, in the case of a generic irrep Rq , 
the mismatch between the patterns ( |5.9| ) and ( |5.16| ) includes additionally the 



perimeter-type contribution owing to the different powers of i^Hi|R,J?o • ^^ 
it is clear from the above discussion, this piece of the answer is irrelevant 
for the determination of the universality class (even for a finite A^ ) which is 
consistent with the previous consideration based on the identification ( 5.15 ). 



6 The large N pattern of the loop-averages. 

Building on the readily computable case of the D = 2 eigenvalue-systems 
( |1.6| ), let us now discuss the general pattern of the large A^ scaling including 
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self-intersecting contours and multi-loop averages. We will argue in partic- 



ular that, in the gauge theories induced from (|1.6| ), the standard large A^ 
factorization takes place. The scaling of a generic irreducible interaction be- 
tween a number of Wilson contours is also considered. In conclusion, we 
comment on peculiar reasons behind the unconventional N —^ oo scaling 



lim < Wl >~ 0(Ar2-^c) ^g^^^ 

of the average of a single nonself-intersecting loop ( |1.7|) . 

It is instructive to begin with the pattern of the irreducible interactions 
between a given number of nonintersecting, nonself-intersecting Wilson loops 
Cp, p = 1, ..., B, on a particular 2d discretized surface of the area L^ . As 
usual, the corresponding amphtudes are defined by the irreducible correlators 
of Bth order 

« n Wi »=< IlWi>-U<Wi> -■■■ ■ (6.2) 

p=i p=i p=i 

They are deduced from the ordinary B -loop correlators < 11^=1 ^c > '^i^ 
subtraction of all lower order reducible contributions as prescribed by the 
standard cluster expansion. One may presume that the pattern ( |4.3| ),( ^^ is 
universal for a wide class of the lattice gauge theories (and in particular for 
the family (|1.1| ) with e~^^^^ given by ( |1.20| )) with the free energy ~ 0{N'^) . 
Then, provided the absence of 'unnatural' cancellation between individual 
diagrams, the 't Hooft double-line representation (for the large A^ WC series) 
yields the topological expansion in the form 

«I[WL»= E w{B,H)N'-'"'\ (6.3) 



Hez 



>0 



where w{B,H) ~ 0{N^) and H denotes the genus of the abstract surface 
in the index-space. The same pattern is supposed to arise in the context of 
the large A^ strong coupling expansion on a lattice. 

To deduce the large A^ pattern of << 11^=1 ^c >> i^ ^^e induced 
gauge theories, I propose to compare the finite N averages < 11^=1 ^c > 
in the ordinary 2d gauge system (|1.1| ) and in the eigenvalue-model (|1.6| ) asso- 
ciated via ( |5.15| ). To this aim, one is to employ the cutting-gluing technique 
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0, H already used in the previous subsection. For simplicity, we assume ad- 
ditionally that (nonself-intersecting) Cp are not homotopic to the defining 
uncontractible cycles of the surface on which both of the associated sytems 
are defined. In this case, the conventional pattern ( |6.3|) gets modified (as it 
will be clear in the continuum limit) by the perimeter-type factors 

«nW^4»=(n^'"''^'') E w{B,H) N'-'^-^ , (6.4) 

P=i p=i //ez>o 

with w{B,H) ~0(A^°). 

To demonstrate ( |6.4|) , we first cut the surface along the set {Cp} of the 
contours that results in a set of disjoint 2d windows. Let the gth window 
have the area Ag being endowed with certain number Bg of boundary discs 
and that Gg of handles. Similarly to the simplest B = 1 case ( |5.2| ), one 
evaluates the corresponding _Bq-disc amplitude 

z^:iAg\{uiz^^^)})=j:e-'^''^^^ u n x«(^p.(zif^)) . (e.s) 



R p=l Ap) 



k ^'^p 



This expression is to be confronted with the corresponding B -disc amplitude 
in the ordinary 2d gauge theory 



Zc-MgliUcJ) = Y.idimRf-'^^-^" e-^(«)^' \{ xr{Uc,) , (6.6) 

R p=l 

where K{R) is defined by eq. (|5.15| ). 



By the same token as in eq. ( |5.5D , the i?-loop correlator < np=i ^c > 
is to be composed from the associated amplitudes ( p. 51 ) (or ( |S.6D respectively) 
performing the remaining Haar integrations over {[/(zj^^ )} . Apparently, this 
procedure can be visuahzed then as gluing the above windows back so that 
their boundaries sandwitch properly the loops Wq involved. Thus, the 
only difference between the ordinary 2d gauge theory and the associated 
2d eigenvalue-system is in the mismatch between the patterns of the multi- 
disc factors (|6.5|) and (|6.6| ). The most transparent situation arises in the 
continuum limit ( |5.20|) when all ratios Ag/ L?' are kept finite. A minor mod- 



ification of the arguments (see Section 5.1) suggests that, in this limit, the 
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identification ( |5.15| ) ensures the proportionality of tlie finite N averages 

< n Wi >e^,.—^ (11 N'-"^'^) < n Wl^ >,au,e (6.7) 



evaluated in terms of the associated multi-disc amplitudes ( |63|) and 
respectively. Presuming the 'smoothness' of the continuum limit, one may 
expect that (|6.4| ) remains valid for generic values of Aq and the coupling 
constants. In particular it implies that in the eigenvalue-systems, despite the 
unconventional perimeter-dependence (|6.1|) , the large N factorization of the 
loop-correlators remains present. 

Turning to the average associated to a single self-intersecting contour, 
we first consider self-intersections on a finite number of sites. The simplest 
example is the eight-figure loop ('two-leaf flower') C^^^ = Ci(xo) U C2(xo) 
composed of the two nonself-intersecting contours Cfc(xo), k = 1,2, which 
share the single site xp in common. Similarly, one constructs a contour 
(7(™) = U^]^Cfc(xo) with the topology of m-leaf flower based at xq . Em- 
ploying the eigenvalue-representation ( ^.7|) , one observes that for such con- 
tour the loop-average is not an irreducible correlator 

m 
fc=l 

in contrast to the pattern of conventional gauge theories. Upon a reflec- 
tion, the large A^ scaling (|6.8|) holds true for the average of an arbitrary 
loop ( |1.7| ) which after removing all the sites with self-intersections is split- 
ted into m disjoint nonself-intersecting contours Ck ■ Complementary, af- 
ter subtraction of all the higher-order 'reducible' parts (contributing as ~ 
0(A^^'^^^^fe=i^^fe), p > 0) the irreducible correlator ( |6.2[ ) of mth order is 
supposed to obey (|6.1|). In particular, for m = 2 we expect that 



1 ^ 

Ut_iCfc{xo) AT i-^ C'fe(xo) I ^ J \ J 



k 



In the case of a contour with self-intersections along links, the associated 
average generically is not irreducible either. The additional reason here is 
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the existence of the (irreducible) one-hnk correlators (with V-Rfc € ^Uk^ ) 

2 

<XRAUp{'z))XR2{U^{^))>-0(Y{dimRk) , ni=n2modiV, (6.10) 



fe=i 



which depend nontrivially on the relevant coupling constants {6^} . As a 
result, the scaling ( |6.1| ) is generically modified as well. 



D 



6.1 Comments on the local large N [Zn] -symmetry. 



Next, let us return to the interpretation of the unconventional pattern (|6.1|) 
of a nonself-intersecting loop average in the gauge theories induced from 
the eigenvalue-systems ( |1.6| ). From the viewpoint of the generic pattern 
73|), ([4.4|) of the matrix theories with the 0{N'^) free energy, the scaling 
]1]) indicates the additional N'^~^^ -suppression. In the more general case 
[l]); the corresponding irreducible average is down in magnitude by the 
factor (n^=i N^~^'-^p) ■ One might expect that (in the WC phase) the latter 
is the result of a specific cancellation between different Feynman diagrams 
(each scaling individually in compliance with the standard pattern ( |6.3| )). As 
this cancellation is not manifest within the formal WC series, we have to 
return to the specifics of the original eigenvalue-models. As we will see, in the 
framework of the saddle-point (SP) method, the extra (0^=1 A^^"^'^'') -factor 
can be traced back to the otherwise unobservable local [Zn]''^ -symmetry 
( p.ll| ) of the leading 0{N'^) order of the effective action S'e//({co'(°)}) (see 



eq. (U)) restricted to the SP orbit {a;(°)} 

To this aim observe first that, were not the invariance (|2.11| ), in the aver- 



age of the product ( |1.8|) each trace could be substituted in the large A^ limit 
by its SP value x/(co'(°)(pfc|zk)) (modulo ( |1.5|) and the Zn gauge transfor- 
mations ( |1.4| ) both invisible for a closed loop C). Altogether, it would result 
in the conventional O(A^) -scaling and trivial perimeter-law pattern of the 
large A^ loop-average ([1. 71). It is the additional averaging over the [Zn]^- 
orbit ( |2.11|) which ensures that (in the eigenvalue-systems (|1.6| )) the leading 
0{N) contribution to < PV^ > is exactly cancelled. Indeed, assume for a 
moment that one would neglect completely the subleading 0{N^~^), /5 > 0, 



orders of S'e//({cij(0)}) in eq. (|2.6|) . In this 'approximation', the fiuctuations 



(in the background of the SP orbit uj^^\pk\'L\^ generated by (|1.4|), ( p..5| ) and 
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( p.lOp ) induced by the infinitesimal transformations ( p.ll| ) become exactly 



zero modes. The contribution of these modes would render an arbitrary 
correlator ( |1.8|) identically zero for any nonself-intersecting loop C 

(6.11) 



{zkGLc} t{zk)eZ^f 



because such Wq does not contain the singlet component with respect to 
the [Zn]^ -averaging over t(zk) G Zn . 

The inclusion of the suhleading orders of 5'e//({co'((/))}) breaks (|2.11|) 
down to (|1.4|), (|1.5|) which makes the above modes quasizero. In turn, it 



justifies that the large N loop-averages (associated to a nonself-intersecting 
contour) are nonvanishing. On the other hand, the cancellation ( |S.ll ) in- 



deed foreshadows the A^^ ^^' -suppressed scaling (|6.1|) . To make it manifest, 
consider the 2-loop average 



< \Wl;\^ >~ 0{N^) , (6.12) 



with iH^cP being invariant under ( p. 11 ). As it is clear from the SP method 



the leading 0(A^^) -order of (|6.12|) is accumulated (in contrast to the ordi- 
nary gauge theories) by the reducible part composed of the rii = n2 = 1 
correlators (|6.10|) . This substantiates that the N —>■ oo average ( |1 . 7|) is 
necessarily suppressed compared to the conventional B = 1 O(A^) -pattern 

m. 



7 Conclusions. 

In this paper we have proposed to induce a class of SU{N) lattice gauge 
theories from the novel family of the vector-field eigenvalue- mo dels (|1.6| ) (or 
their generalizations) employing the mapping ( |1.9|) . The latter models, in 
addition to the local [U{N)]^ conjugation-invariance ( |1.3| ), are endowed with 
the Zn gauge symmetry ( |1.4|) and global [Zn]-^ invariance (|1.5| ). The crucial 
consequence of ( p..3| ) is that both the action (|1.6| ) and the Wilson loop 
averages (|1.8|) depend nontrivially only on the eigenvalues a;(pfc|zk) of the 
relevant link- variables. 
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Consequently, both the partitition function X^d (composed of ( |2.1| )) 
and the loop-averages (|1.8| ) can be rewritten in terms of the GLR coeffi- 
cients ( |2.2| ) parametrized by 0{N) integers. In turn, the 1/N expansion 
of both the free energy and (|1.8| ) can be reformulated as the semiclassical 
series provided that —InlXx^o] ~ 0{N'^) . In this perspective, the eigenvalue- 
models ( |1.6| ) facilitate a master-field representation of the associated induced 
gauge systems. Recall also that in any literally continuum system the local 
conjugation-invariance ( |1.3| ) is not consistent (even at the classical level) 
with the presence of any finite number of derivative- dependent terms in 
the action. Thus, from the viewpoint of the continuum YM theory, the 
invariance ( |1.3|) of (|1.6|) is to be interpreted as the huge auxiliary symme- 
try introduced by the judicious discretization (regularization) of both the 
action and the space-time. 

In the large A^ limit, the partition function of (|1.6| ) is reproduced (via 
( |1.14| )) by the (L^th power of the) reduced generating functional ( p..l7| ) of 
the GLR coefficients which can be further transformed into the 1-matrix 
representation ( p..l9| ). To be even more specific, the 1-matrix eigenvalue- 



model in ( |1.19|) is selected in the simplest form ( |1.20| ) which makes accessible 
the complete large N phase structure of the D >2 induced gauge theories. 
In particular, we pay the special attention to the issue of the continuum limit 
(CL) in the latter theories. The important place in this analysis is played by 
the scaling-condition ( |4.6| ), imposed on the effective 1-matrix system ( |1.19|) . 
As it is demonstrated, this condition predetermines that in the induced gauge 
system the localization {Up{z) — > 1} of the link- variables takes place which is 
tantamount to CL. Also the phenomenon of the auxiliary 'continuum' limits 
is discussed. 

The computation of the large A^ Wilson loop averages < VF^ > is par- 
ticularly simple in the D = 2 case where the concise representation ( ^.11 



IS 



derived for generic nonself-intersecting loop C on an arbitrary 2d surface. 
Building on the transparent relation with < Wq > evaluated in the ordinary 
2d gauge theories, we reveal the peculiar modification ( |6.4| ) of the standard 
pattern ( |6.3| ) of the 1/A^ -topological expansion. However the physical mass- 
scale (like the string tension in (|5.11| )) is adjusted to scale as 0{N^) , and the 
large A^ factorization is expected to remain present. Therefore, we conclude 
that the additional N^'^c: -suppression ( p.l|) does not prevent to consider 
the induced gauge systems as the elligible large A^ cousins of the conven- 
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tional SU{3) lattice gauge theories. Moreover, it is reasonable to expect 
that the 0{N^) physical string tension (if any) in the induced lattice YM 
system reproduces that in the associated continuum large N YM theory. 
It is supported, in particular, by the explicit D = 2 derivations of Section 
5.1. 

Definitely, the most ambitious goal of our project is to find an approach 
for computation of the D = A large A^ loop-averages in the proposed in- 
duced systems adjusted to have the proper continuum limit. If successful, it 
might provide with the framework appropriate to address the issue of con- 
finement in the standard D = 4 continuum SU{N) theory in the limit 
N ^ oo . Unfortunately, the D = 2 technique of Section 5, being directly 
generalized to D > 3 , does not seem to provide with a practical scheme. 
More promising direction is to synthesize the latter technique with the 1/A^ 
saddle-point method applied to the eigenvalues u!{pk\z]^) employing the ap- 
propriate representation ( |1.8| ) of liniTv^oo < M^c ^ • '^^^ actual analysis in 
this framework is additionally complicated by the presence of the quasizero 
modes associated to the [Zn]^ symmetry ( |2.11| ) of the leading 0{N'^) order 
of the effective action S'e//({ct;'^°^}) (defined by eq. ( p.6|) ) restricted to the 
SP orbit {u;*^°)} . As a result, even in the computation of limjv_^oo < W^c > 
it looks necessary to take into account the suhleading orders of Seff{{oj}) . 
Nevertheless, one might expect that the presumable large A^ factorization 
( |6.4|) foreshadows potential simphfications. Clearly, it calls for a new piece 
of technology to be developed. 
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Appendix A: The local [Z^v]^ - invar iance. 

In this appendix we discuss the local [Zjv]^ -symmetry ( p.ll| ) of the leading 
0{N'^) order of the effective action Seff{{uj^^^}) (defined in eq. (p^)) re- 
stricted to the SP values {uj^^^ . Owing to S{D{D — l)/2) invariance ( |1.13|) , 
it suffies to consider the simpler version of the SU{N) plaquette-factor ( |1.6| ) 
with the sum running over the single species R2 

z{{ujm)=j:e^'^''^^ n xr.{u^{m)-xrM'^\^))^ (a.i) 

R2 {t^^} 

where S{R2) = S{R2) and (p E 2® p . 

Our purpose is to demonstrate that the SU{N) SP solution uj^'^\(j)) is 
supposed to be unique and Z2-selfdual (where xni.'^i.^')) = [xr{.^{4>))]^) 

u;f\4>) = -ooSl.^M ^ cuf (0) mod S{N) ® [Z^]^ , (A.2) 



modulo the Weyl ^(A^) group ( |2.1CI| ) combined with the [Zn]-^ transfor- 



mations ( p.lll) . To this aim, one first observes that the factor (|A.1|) is by 
construction Z2 self dual Z({u;(0)}) = Z~^{{uj{(j))}) = Z{{uj{(j))}) . It implies 
that on the SP orbit, generated by the relevant symmetries to be determined, 
there is a point corresponding to a Z2 self dual (modulo ( p.lO|) ) solution 
{a;'^°)(0)} . Consider the SP stt6orbit generated from the above Z2 invariant 
solution by the (local) Weyl group ( p.lO| ). In this case, the corresponding 



stronger version of the constraint ( [A .21) is tantamount to the reality of the 
character XR(a;(0)) = [xR{u){(f)))]'^ for any (not necessarily selfdual) irrep 
R . In turn, the latter property is sufficient for the selfconsistency of the local- 
ization ( |3.1D on R2 = R2 ■ otherwise the effective action 5'e//(-R2|{^(0)}) 
for R2 

^e-5.,,(^.,M,)})^^(^^(^)j) (A.3) 

R2 

would not be Z2 -selfdual with respect to the interchange R2 ■^^ R2 ■ 

Now we are in a position to prove [Zn]'^ transformations ( |2.11| ) of the 
above S{N) suborbit {uj^'^\(J))}s(n) ■ Let us introduce the twisted action 



Seff{R2, t\{u;{(j))}) , resulting (akin to ( |A.3D ) from eq. ( [A .11 ) after the [Zn 



^D 



deformation of the characters Xr(^) involved 

?i ^ '7 _ D ^ \ 

■n(i?) 

42 



Xr{V) ^ XRitV) = e'^^'hRiV) ; t = e''^eZ,,, Re F^, (A.4) 



(where n{R) is defined by (^) so that Z{{uj{(p)}) -^ Z({cj(0)}, t) . We 
claim that Seff{R2,t\{^^'^K<P)}) remains to be Z2-selfdual 



5e//(i?2,t|{u;W(0)}) = %;(i?2,t|{c.W(0)}) (A.5) 

provided that {cu^^^cj))} E {lu^^\(P)}s(n) , i-e. fulfils the stronger version 
of ( |A.2| ) where the [Zn]'^ group is omitted. Moreover, we assert that (for 
such {ci;(°^(0)}) the resulting from (|A.4| ) sum (|A.3| ) over SU{N) irreps R2 
is localized for Vt G Zn on the same i?2 (t) = R2 (1) = R2 as before 
[Zn]-^ -twisting. In turn, it is evidently tantamount to the invariance of the 
leading 0{N'^) order of ^e//({^^°^}) (defined by eq. (p|)) under (|ll|). 



To prove these statements, one notes first that the explicit form ( p.7| ) 
of Z2 transformation yields exp{n{R)ln[t]) = exp{—n{R)ln[t]) . As only 
the real part of the overall twist-factor contributes into Z{{uj{(f))},t) , the 
latter identity justifies ( |A.5D . In the SU{N) case (where A^ =0 



see (^4.22|) ) it implies that the saddle-point irrep is supposed to be selfd- 



ual 4°^(t) = 4°^(t) irrespectively of t . To justify i?f (t) = i?5°^(l), we 
employ the elligibility to restrict the large A^ twisted sum (|A.3|) to that 
over SU{N) selfdual irreps R2 = R2 for which A^Ai = 2EiIiAi G [2Z>o] 
so that ^'^(^a) — exp[i7i{Xi — A^ + 1)] . Finally, combination of the latter 
equation with the choice A^ G [2Z>o + 1] results in Ai G [2Z>o] so that 
^n{R2) _ -^ Presuming the 'smoothness' of the large A^ limit, it substanti- 
ates t -independence of R2 (t) . 



Appendix B: The D = 2 loop- averages. 



Consider the large N D = 2 average ( |5.9| ) of a nonself-intersecting Wil- 
son loop in the case when Rq = f and the weight e"'^'^^-' belongs to the 



U{N) subvariety ( 4.26 ) of the pattern ( 1.20 ). To begin with, as we will see 
in a moment the factor BL^{R\f,A) scales as 0{N'^~'^'^) provided the in- 
tegers \j (defining i?({A}) ) are of order of ~ 0{N) in the limit N ^ oo . 
We assert that the large A^ sum in eq. ( ^.9|) is localized on the solution 
i?({A(°)}) = R{X^^\t)) of the SP equations (j^). Indeed, recall first that 



a (unique, Z2 -selfdual) R^'^^ = R{{\^^^}) is the irrep on which the large A^ 
sum ( ^.41) (defining the partitition function) is localized. It is reasonable to 
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presume that, for |A^- — Xj\ ~ 0{y/N) , there is such a > that 



i?L,(i?({AW})|/,A) 



OiN-'^) , (B.l) 



i.e. BL^{R{{\})\f,A) behaves smoothly on the scale of the characteristic 
fluctuations of Xj . Altogether, it justifies the above assertion. 

Next, one is to derive an explicit form of Bi^{R^^^\f^ A) in terms of the 
associated spectral density ( |4.25| ). The representation theory p tells that 
L^^ 1^(0) J is nonzero when the associated to Ri (double) Young table can 
be obtained from that corresponding to -R*^°^ by adding a single box in an 
admissible way. The latter addition is elligble into the i th row if and only 
if the number of boxes rij in this row (see eq. (|3.5| )) is strictly less than in 
the preceeding one 

rii < n,_i ^^ (Ai_i - A,) > 2 . (B.2) 



According to the definitions (|4.21|) and (|4.25|) of dX{t)/dt and p{ri) , it implies 



that in the domain of the coupling constant(s) where < p{rf) < 1/2, Vr], the 
box can be generically added to any row of the Young table corresponding 

to pirj) = {dt/dh^^\t))\^{o)(^'^^^ . Indeed, owing to ( |B.1D one can choose any 



discretization Xj /N, j = 1,...,N, of the continuous function X^^\t) = 

/i(o)(t) + l/2 which satisfies |A^A(o)(j7A^) - Af ^| « VN . Condition p{r]) < 
1/2 evidently ensures the existence of at least one such discretization obeying 
T2I1 



Due to ([4.30|) , the constraint p(?7) < 1/2 is always fulfilled in the large 
N WC limit {bj. — > 0} relevant for the analysis of the continuum limit. 
Therefore in what follows, we will restrict our attention to the cases when 
piv) ^ 1/2 and thus ([B.2|) ) is valid. As a result, (by the same token as in 
|T^) one obtains 

Ee-'^^^^)-^^^"'^l^^K,^(o),/ = [dv Piv) [nv)r e-'^\ (B.3) 

k ^-1 1 /N 

T{h{-)) = expij: ln{l + ,/ )] , (B.4) 
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where we have skipped the upperscript: h^^\t) = h(t) . In the derivation of 
the above equation, we employ that when the box is added in kth row (of 
the Young table associated to R^'^^ ) then 



Mo rx(0) _ N^i2m-l ^V(n] 

m=l 



N2m-1 ^^ lr?=M#)' 






where V{R{{X})) is defined by eq. (|;^. 



The remaining step is to find the continuum representation for T(?7) in 
terms of p{ri) . For this purpose, one is to expand the logarithm in eq. ( p.4|) 

nxj^)) - -( E + E ) E -[iv(/4) - M^))l- . (B.7) 

if^k \i-k\>M p>l -f^ 

treating {h{-^) — h(j^)) differently for small, \i — k\ < M , and large, \i — k\ > 
M , values of \i — k\ (where M = N'^, < 7 < 1 ). In the former domain, 
one is to substitute 

while in the latter case the (partial) sum over i is to be transformed into 



As it is clear from (p.9| ), the situation is somewhat distinct for the p = 1 
and the p > 2 terms of the expansion ( [B.7| ). In the p > 2 case, the 'long- 
wavelength' term ( p.9|) does not contribute into the leading 0{N^) order 
(of ln[T{ri)] ) we are interested in for eq. (|B.3|) . According to ( p.8|) , the 
remaining 'short-wavelength' p > 2 terms are represented by the convergent 
sums which results (see [|l4l) in 
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where h{^) = r] and the contribution of odd p = 2m + 1 vanishes. Com- 
plementary, in the p = 1 case the 'short-wavelength' contribution goes to 
zero while the large A^ limit of the 'long- wavelength' part merges with the 
principle value of the integral P J _f . Combining all the pieces together 
and employing the SP equations ( |4.27|) , one arrives at ( p.ll|) . 
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